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1 Introduction and statement of results 

In a large number of problems coming from statistical or quantum mechanics, 
involving real partial differential operators of order two with the spectrum 
contained in the right half-plane, one is often interested in the splitting be- 
tween the two smallest real parts of the eigenvalues, at least when it is possi- 
ble to show that the first one is simple and isolated. This kind of issue may 
concern the Schrodinger operator or the Witten Laplacian in quantum me- 
chanics, and the Kramers- Fokker- Planck operator, or more generally, some 
models of chains of oscillators coupled to heat baths, in statistical mechanics. 

In the semiclassical context, the tools available for studying the splitting are 
very powerful, and allow for a detailed analysis for models where a natural 
small parameter is in play, e.g. the Planck constant or the low temperature 
parameter. In some cases the splitting may be exponentially small with 
respect to the parameter, and a so-called tunneling effect may appear, related 
to a very low rate of convergence for the associated evolution problem. It is 
often related to some degenerate geometry in the model, such as the presence 
of multiple confining wells. 

In fact, for some of the preceding semiclassical differential operators, a par- 
ticularly convenient structure is sometimes available, which simplifies consid- 
erably the analysis of the splitting and of the tunneling effect. Indeed, when 
the operators have a Hodge Laplacian type structure of the form P = d*d, 
where d is a (perhaps modified) de Rham operator, the eigenspaces have 
some natural stability properties with respect to the operator d, and the 
study of the eigenvalues can be reduced to the study of the singular values 
of d on finite dimensional subspaces. This type of operators is sometimes 
called supersymmetric (see in particular [2TJ and the full definition given be- 
low), and this method was successfully employed for determining the splitting 
corresponding to some models mentioned above. 
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The natural question raised in this article is to determine under which condi- 
tions a given (not necessarily selfadjoint) operator is of supersymmetric type. 
The answer given here implies that for some models describing systems out 
of equilibrium, there is no such property, and we give a concrete example in 
the last section of the paper. The significance of this kind of result is that 
the study of the splitting for these models may be far more complicated, in 
particular in the non-self adjoint case, and to the best of our knowledge, no 
such analysis has been performed in the non- self adjoint case so far. 

By way of introduction and also, for motivation, let us mention now some 
examples of supersymmetric differential operators and related results about 
the eigenvalue splitting. The first example is the Witten-Hodge Laplacian, 

W = -h 2 A + (d x V(x)) 2 -hAV(x), (1.1) 

on R™. Here V e C°°(R n ;R) is a Morse function with d a V e L°°(R n ) for 
|a| > 2, and such that |9 x ^(x)| > 1/C when \x\ > C > 0. The second 
example is the Kramers- Fokker- Planck operator given by 

K = -h 2 A y + y 2 -hn + y-hd x - d x V(x) ■ hd y , (1.2) 

on R^ y , where the potential V has the same properties as above. This is a 
kinetic (non-selfadjoint) model for an oscillator coupled to a heat bath. 

For these two models, the splitting can be evaluated very precisely in the 
semiclassical limit h — > 0. Considering W and K as unbounded operators 
on L 2 , it was shown in [9], [8] for the Witten Laplacian, and in [10], [11] 
for the Kramers-Fokker-Planck operator, that, under suitable assumptions, 
it is exponentially small. We notice here that in both cases, the operators 
have non-negative symmetric parts and have \L\ = as a simple eigenvalue, 
under the additional assumption that V(x) — > oo as x — » oo. For simplicity, 
let us recall here the result in the case when the potential V has two local 
minima m± and a saddle point sq of signature (n — 1, 1), although a far 
more general result is given in [12] . in the Kramers-Fokker-Planck case. Let 
S± = V(sq) — V(m±) and S = min(5 + , SL). Then there exists c > and 
ho > such that for all < h < h there are precisely two eigenvalues of K 
in the open disc D(0,ch), one of them being /ii = 0, and the second one is 
real positive, given by 

fi2 = hl(h)e~ s/h , (1.3) 
where 1(h) = l + hl\ + h 2 l 2 + h 3 l 3 ..., with l > 0. 

We shall now discuss the notion of a supersymmetric structure more precisely. 
Let X be R™ or a smooth compact manifold of dimension n, equipped with 
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a smooth strictly positive volume density co{dx). We can then choose local 
coordinates x±, .., x n near any given point on X such that co{dx) = dx\-...-dx n . 



Let d : C°°{X; A k T*X) -> A fc+1 T*X) be the de Rham complex, and 

let 5 : C°°(X; A k+l TX) -> C°°(X; A k TX) be the adjoint of d with respect to 
the natural pointwise duality between A k T*X and A k TX, integrated against 
co. In the special local coordinates above, we have 



Let A(x) : T*X — > T X X be a linear map, depending smoothly on x G X . 
The /c-fold exterior product A k A maps A k T*X — > A k T x X, and by abuse of 
notation, we shall sometimes write simply A instead of A k A. By convention, 
A A is the identity map on R. Associated with A, we have the bilinear 
product on the space of compactly supported smooth fc-forms, 



{u\v) A = (Au\v) = / (A(x)u(x)\v(x))co(dx)i u,v G C^(X; A k T*X). (1.5) 



Since we shall restrict the attention to real-valued sections and operators, 
there is no need for complex conjugations in (II. 5ft . When the map A{x) is 
bijective for each x, there is a natural way of defining the formal adjoint of 
a linear operator B taking fc-forms to £-forms on X. This adjoint, denoted 
by B A '*, maps £-forms to fc-forms and is given by 



Here B* is the usual adjoint taking ^-vectors to /c-vectors. Notice that when 
k = 0, we can still define the adjoint B A '*, even when A is not everywhere 
invertible. In the case when B = d, the de Rham differentiation, we get 




(1.4) 



(Bu\v) A = (u\B A '*v) A , 



and more explicitly, by the following expression, 



B A '* = (A k A t )- 1 B* A £ A\ 



d A >* = (A'y^A 1 , 

and for the adjoint of the restriction of d to 0-forms, we get 

d A >* = SA\ 



We can write, in the special local coordinates above, 
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Let ip G C°°(X;R) and let us introduce the Witten complex, given by the 
weighted semiclassical de Rham differentiation 

d vh = e -^ h ohdoe v/h = hd+{d<p) A : C °°(X; A k T*X) -»• C °°(X; A fc+1 T*X). 

(1.6) 

Considering the adjoint of d^h with respect to the bilinear product (ll.Sp . we 
can introduce the operator, 

d%* h = e <p/h ohd A '*o e - <p/h , (1.7) 

which on 1-forms becomes, 

d^ = (hS + d^)oA\ 



We are now able to give a precise definition of the notion of a supersym- 
metric differential operator. In doing so, it will be convenient to distinguish 
between the non-semiclassical case, corresponding to the situation where the 
parameter h > is kept fixed, say, h = 1, and the semiclassical case, where 
we shall let h G (0, ho], ho > 0. 

Definition 1.1 (i) (The non-semiclassical case.) Let P = P{x,D x ) be a 
second order scalar real differential operator on X, with C°° -coefficients. 
We say that P has a supersymmetric structure on X if there exists a linear 
map A(x) : T*X — > T X X , smooth in x E X , and functions ip,ip G C°°(X; R) 
such that 

P = d£*d v . 
Here d^ = and similarly for d^. 

(ii) (The semiclassical case.) Let P = P(x, hD x ; h) be a second order scalar 
real semiclassical differential operator on X . We say that P has a super- 
symmetric structure on X in the semiclassical sense if there exists a linear 
h-dependent map A(x; h) : T*X — > T X X , smooth in x G X , and functions 
ip = tp(x;h) and ip = ip(x;h), smooth in x, such that ip = tpo(x) + O(h), 
if) = ipo(x) + 0(h) in the C 00 -sense, and for which we have 

P = d$; h d v>h , (1.8) 

for all h G (0, h ], h > 0. 

Remark. Let us notice that no control on the behavior of A(x; h) is required, 
as h — > 0, in the semiclassical case in Definition 1.1. It is an interesting 
problem to determine when we can find A(x; h) with some uniform control 
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as h — > 0. As we shall see from the proof of Theorem 11.21 the symmetric 
part of A is immediate to determine, while the antisymmetric part arises as 
a solution of a 5-problem with exponential weights. 

Assuming that the linear map A(x; h) is invertible for each x G X and all 
h G (0, ho], let us put in an arbitrary degree k > 0, 

U A,il>,<p = d i,',h d <p,h + d^hd^lh- (!- 9 ) 

The supersymmetric semiclassical operator P in (II. 8p agrees with the re- 
striction of ^A,ip,ip to 0-forms. The analysis of the eigenvalue splitting of P 
strongly depends on the following formal intertwining relations, 

n ( fc+1 ) / 7 —A n( fc) jA,* n (k+l) _ {k) jA,* /i in\ 

LJ J 4,Vw u< M ~~ u <M LJ J 4,i/w u iM A,f,(p ~ A,tl),<p u tp,h- l ± - ±u J 

As a clarifying example, let us now come back to the case of the Witten 
Laplacian W in (II. ip . A straightforward computation shows that the op- 
erator IV enjoys a semiclassical supersymmetric structure in the sense of 
Definition 1.1, 

W = {hd+ {dV) A Y'*(hd + {dV) A ) = d\* h d Vth , 

with tp — if; — V, and A — I. In the case of the Kramers-Fokker-Planck oper- 
ator K, the semiclassical supersymmetric structure was observed in pTJ, [2], 
and in [20]. It will be recalled in the next section, following [TO] . 

The first result of this paper is a necessary and sufficient condition for a 
second order real differential operator to be supersymmetric, both in the 
non-semiclassical and in the semiclassical sense, either locally or globally on 
X. 

Theorem 1.2 (i) (The non-semiclassical case). Let P be a second order 
scalar real differential operator on X , which can be written in local coordinates 
as follows, 

n n 

P = - ^2 dx i ° B iA x ) od x k + } y Vj(x) o d X] + v , 

j,k=l j=l 

where (Bj^) is symmetric, and (Bj ^), Vj, Vq are real-valued and smooth. In 
order for P to have a supersymmetric structure on X it is necessary that 
there should exist ip,ib G C°°(X; R) such that 

P(e"^)=0, P*( e -^)=0. (1.11) 
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Here P* is the formal adjoint of P with respect to the L 2 -scalar product deter- 
mined by the density of integration u(dx). If (11. lip holds and the 5-complex 
is exact in degree 1 for smooth sections, then P has a super symmetric struc- 
ture in the sense of Definition 1.1. 

(ii) (The semiclassical case.) Let P be a second order scalar real semiclas- 
sical differential operator on X , which can be written in local coordinates as 
follows, 

n n 

P = - ^ hd Xj o B j!k (x; h) o hd Xk + ^ Vj(x; h) o hd Xj + v (x; h), (1.12) 

j,k=l j=l 

where {Bj^) is symmetric, and {Bjk), Vj, v are real-valued, h-dependent, 
smooth in x G X . In order for P to have a supersymmetric structure on 
X in the semiclassical sense it is necessary that there should exist functions 
<f(x; h),ifj(x; h) smooth in x G X, with ip = ipo(x) + 0{h), if) = ipo(x) + 0{h) 
in the C 00 -sense, such that 

P( e -^ h ) = 0, P*(e-^ /h ) = 0, (1.13) 

for all h G (0, ho], h > 0. If f 1 1 . 1 3 j) holds and the 5-complex is exact in degree 
1 for smooth sections, then P has a semiclassical supersymmetric structure 
in the sense of Definition 1.1. 

Remark. Assume that X is a compact Riemannian manifold and that P = 
— A + v, where v is a smooth real vector field. Then (11.111) holds, with 
tp — if; — 0, precisely when the vector field v is divergence free. 

As a fundamental example of a second order semiclassical operator for which 
the supersymmetric structure may break down, we shall now discuss the case 
of an operator associated to a model of a chain of anharmonic oscillators 
coupled to two heat baths. Referring to the next section for the construction 
of this model, let us consider here the following real semiclassical differential 
operator of order two, 




+ yhd x -{d x W(x)+x-z)-hd y . (1.14) 

Here x = (xi,x 2 ), y = (2/1 , 2/2) ? and z = (zi,z 2 ) belong to R 2n . Furthermore, 
7 > 0, atj > 0, j = 1, 2. The effective potential W is of the form W(x) = 
W 1 {x 1 ) + W 2 {x2) + SW{x 1 ,x 2 ), where Wj G C°°(R n ;R), j = 1,2, and 5W G 
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C°°(R 2n ; R). The parameters ola are proportional to the temperatures in the 
baths, j = 1,2. As we shall see, the structure of operator Pw is completely 
different depending on whether «i = «2 or not. The following is the second 
main result of this work — see also Theorem 13.11 below for a more precise 
statement. 

Theorem 1.3 Consider the operator Pw defined in fll,14p . We have the 
following two cases. 

i) If oi\ = «2 (the equilibrium case) or 5W = (the decoupled case), then 
Pw has a semiclassical supersymmetric structure on R 6n , in the sense 
of Definition 1.1. 

ii) Take 7 = 1 and assume that ct\ 7^ cti- Then there exist a Morse function 
W\[x\) with two local minima and a saddle point, a positive definite 
quadratic form ^2(2:2), and SW(xi, £2) £ C^(H 2n ), arbitrarily small in 
the uniform norm, such that the operator Pw has no supersymmetric 
structure in the semiclassical sense, on R 6n , in the sense of Definition 
1.1. 

The study of the tunneling effect for the Witten Laplacian through its super- 
symmetric properties was first performed in [9J, following the supersymmetry 
observation pointed out in [21]. As already mentioned, this structure is of 
great help in the study of the eigenvalue splitting, since in particular it al- 
lows one to avoid the study near the so-called "non-resonant wells", such as 
the saddle point in the example mentioned above. The results were subse- 
quently generalized in [3], [I], in [8] with a full asymptotic expansion, in [II] . 
[15] in cases with boundary, and in [16], [IT] in the case of forms of higher 
degree. Notice that the computation of the exponentially small eigenvalues 
is performed here using the singular values of the Witten differential. 

For the Kramers- Fokker- Planck operator, the supersymmetric structure was 
first observed in [TJ, [2], and [2D], and was used for a complete study of the 
tunneling effect in [10] and [12]. This structure helps substantially in this 
non-self adjoint setting, and using also an additional symmetry of a PT-type, 
we were able to compute the singular values, in order to get results in the 
form given in (II. 3p . Let us also remark that if the supersymmetric method is 
not available, a natural approach to the tunneling analysis would proceed by 
means of exponentially weighted estimates of Agmon type, where the notion 
of the Agmon distance should be replaced by a degenerate Finsler distance, 
see PS]. 
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In fact in [TU], [TT], [12], it was shown, under additional assumptions, that 
the tunneling of any supersymmetric operator of order two could be studied 
using the method given there, and the case of chains of oscillators coupled 
to two heat baths was given there as an example. This model was first 
introduced in [6] and was shown to be supersymmetric in the case of the 
equal temperatures in [TT] . For general second order non-negative operators, 
the question of intrinsic geometric structures and their relations with the 
spectral gap in asymptotic regimes was also studied in [7] in the elliptic case, 
and this article was one of the motivations for our study of supersymmetric 
hidden structures here. 

The paper is organized as follows. In Section 2, we first give a necessary and 
sometimes sufficient condition for a general real partial differential operator 
of order two to be supersymmetric, both in the non-semiclassical and the 
semiclassical cases, and establish Theorem 11.21 We then review the three 
supersymmetric examples and in particular, discuss the equilibrium and the 
uncoupled cases for the chains (case i) of Theorem 11.31 Section 3 is devoted 
to the study of the case of different temperature for the chains and to the 
proof of part ii) of Theorem 11.31 

2 Generalities on supersymmetric structures 
and some examples 

In the beginning of this section we shall establish Theorem 11.21 Let us begin 
with the non-semiclassical case h = 1. Let X be R n or a smooth compact 
manifold of dimension n, equipped with a smooth strictly positive volume 
density u(dx). Assume that local coordinates X\, ..,x n have been chosen so 
that u(dx) = dx\ ■ ... • dx n . If 

is a smooth vector field on X, then the divergence of v is well defined by the 
choice of u and in the special coordinates above, we have 

n 

divu = (2.2) 

i=i 

Let d be the de Rham exterior differentiation and let 8 be its adjoint, defined 
in the introduction. Then using the local expression (11.41) . we obtain that 

Sv = — div (v). 
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Let now A(x) : T*X — > T X X depend smoothly on x G X, and let us consider 
the bilinear product defined in fll.5p and the induced action of A on /c-forms. 
We now let P be a second order real differential operator on X, which we 
can write in the special local coordinates above as follows, 

n n 

P = - 5^ d Xj o B jtk (x) o d Xk + ^ v j( x ) ° + u o- (2-3) 

j,k=l j=l 

Here (-P/,fc) is symmetric and Bj^, Vj and v are real-valued and smooth. 
Viewing P as acting on 0-forms, we first ask whether there exists a smooth 
map A(x) as above, such that 

p = d A '*d = 6A% (2.4) 

either locally or globally on X. 

Proposition 2.1 In order to have (12.41) . it is necessary that 

P(l) = and P*(l) = 0. (2.5) 

i/ere P* i/ie adjoint of our scalar operator P with respect to the L 2 -scalar 
product determined by the density of integration u(dx). 

Proof. If U2I1I) holds, then clearly P(l) = 0, since d(l) = 0. Since P* = 5 Ad, 
we also have P*(l) =0. □ 

Proposition 2.2 The property 02 .5p zs equivalent to the following property, 

vq = and divv = 0. (2-6) 
Proof. In the special local coordinates, we see that 

n n 

P* = d Xj o B jjk (x) o d Xk - ^ v i ° 9 x s ~ div v + v . 

j,k=l j=l 

Thus the property P(l) = is equivalent to v = and the property P*(l) = 
is equivalent to Vo — div (v) = 0. □ 

We now look for a smooth map A = A(x) such that (12.41) holds. In the 
special coordinates above, let A = (Aj^(x)) be the matrix of A. Then, 

n 

5A'd = -J2 d *i ° A k,j( x ) ° d * k - 
j,k=i 
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Write A = B + C, where B is symmetric, B t = B, and C is antisymmetric, 
C* = — C. Then we get 

n n n 

»\'<i = - E ^ ° ^•^•'•) ° ^ + ECCm^.*))^- (2.7) 

j',fc=l k=l j=l 

In order to have (I2.4p . we see that -B has to be equal to 5, and we assume 
that from now on, so that 

A = B + C, B l = B, C* = -C, (2.8) 

where 5 is the matrix appearing in (12.31) and the antisymmetric matrix C 
remains to be determined. In order to have (12.41) it is necessary and sufficient 
that 

n 

v k = Y,d Xj (C j)k ). (2.9) 

Consider the 2-vector form 

n 

j,k=i 

By a straightforward calculation, see f ll.4p . we obtain that 




so (12. 9p amounts to 

i; = -lar. (2.io) 

As we have seen, the necessary condition (12.51) is equivalent to (12. 6p . which 
contains the assumption that divv = 0, i.e. that Sv = 0. Since the 5-complex 
is locally exact in degree 1, we get the following result. 

Proposition 2.3 If ( 12. 5 j) holds, then locally, we can find a smooth matrix 
A such that ( 12. 4 p holds. More precisely, we can find a smooth matrix A so 
that (12. 4p holds, in any open subset X C X , where the 5-complex is exact in 
degree 1 for smooth sections. 
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Remark 2.4 If A = B + C is a solution to (I2.4p . then A = B + C is another 
solution if and only if C — C (identified with a 2 -vector form) is 5-closed. 

We shall now replace the assumption ( 12. 5 p by the more general assumption 
that there exist smooth strictly positive functions e~ v and e~^ in the kernels 
of P and P* respectively, 

P(e^) = 0, P*(e~^) = 0. (2.11) 

Put P = e"*oPo e-v, so that P* = e~ v o P* o e~^. Then P satisfies fl2~5]) . 
Hence if X C X is an open connected subset where the 5-complex is exact 
in degree 1, we have a smooth matrix A on X such that 

P = 5A l d in X. (2.12) 

Putting d v = e~ ip odo e v , <fy = e - ^ o d o e^, we get with A = e^ + <M: 

p = e ^P> = e^&T^VW = dJA*tZ v = (c^) A *<V 

Summarizing the discussion so far, we see that we have now established the 
first part of Theorem 11.21 addressing the non-semiclassical case. Restoring 
now the semi classical parameter h G (0, ho] and recalling the notion of a 
supersymmetric structure in the semi classical sense given in Definition 1.1, we 
see that the arguments given above can be applied to the conjugated operator 
e -VM o P(x, hD x ; h) o e~^ h , where up = (p (x) + 0{h), ip = Mx) + 0(h), for 
each fixed value of h G (0, ho]- The second statement of Theorem [L2] follows 
and the proof of Theorem 11.21 is therefore complete. 

Corollary 2.5 Let P = P(x, hD x ; h) be a second order real semiclassical 
differential operator on X, which in local coordinates, can be written as fol- 
lows, 

n n 

P = - 7^ hd x . o B j>k (x) o hd Xk + ^ vj(x) o /i<9 Xj + v (x) + hw (x). 

j,k=l j=l 

Here (Bj t k) is symmetric, and (Bj t k), Vj, vo, wo are real-valued, smooth, 
and independent of h. Then a necessary condition for the operator P to be 
supersymmetric in the semiclassical sense is that there should exist smooth 
h-independent functions <po and ip satisfying the following eikonal equations 
on X , 

n n 

^ B jtk (x)d Xj ipod Xk ipo + y^ y Vj(x)d Xi ipo(x) - v (x) = 0, 
j,k=i j=i 
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and 



n n 

- ^2 B j,k( x )9 Xj ^odx k ipo + ) y Vj(x)d Xi i/;o(x) + v (x) = 0. 

j,k=l j=l 



Remark. Let us notice that if an operator P possesses a supersymmetric 
structure in the semiclassical sense, then so does e^' h oPoe~^ h , for a smooth 
/i-independent function /. Indeed, if P = d^d^h, then 

pf/ h o P o p- f l h - d A '* d t x. 

We shall end this section by coming back to the three examples mentioned in 
the introduction, namely, the Witten Laplacian, the Kramers-Fokker-Planck 
operator, and the model of chains of oscillators coupled to heat baths. All 
three of them come from stochastic differential equations, and we refer to Sec- 
tion 5 in [11] for a complete discussion concerning their derivation. Here we 
shall merely quote some results given there. The associated time-dependent 
equation, 

{hdt + P)f{t,x) = 0, t>0, (2.13) 

describes the evolution of the particle density, and due to the stochastic 
origins of this model, the equation P*(l) = is always satisfied, so that for 
= 0, the equation P*( e -^ /h ) = is given for free. According to Theorem 
11.21 the only remaining thing to check in order to obtain the existence of 
a semiclassical supersymmetric structure for P is then the existence of a 
smooth function cp(x;h) = tfo(x) + 0(h), such that P(e~ 2lfi ^ h ) = 0, for all 
h > small enough. 

The Witten Laplacian. We begin with the semiclassical Witten case. It 
corresponds to an evolution equation with a gradient field —^f'VV(x) and a 
diffusion force coming from a heat bath at a temperature T = h/2. The 
stochastic differential equation corresponding to this model is of the form, 

dx = —•yd x Vdt + y^yhdw. 

Here x 6 R™ is the spatial variable, the parameter 7 > is a friction coeffi- 
cient, and w is an n-dimensional Wiener process of mean and variance 1. 
Equation (I2.13P for the particle density in this case is then 

hdj - ^hd x ■ (hd x + 2d x V)f = 0. (2.14) 

We see that P*(l) = 0, and also that P(e~ 2lp / h ) = 0, for all h > 0. Here 
we have written tp — V. An application of Theorem 11.21 gives that P has 
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a supersymmetric structure on R n in the semiclassical sense. The function 
M(x) = e~ 2 ^/ h is sometimes called the (non- normalized) Maxwellian of the 
system and corresponds to a stationary density, provided that it is integrable 
on R n . 

Writing / = M y l 2 u 1 we obtain from f)2.14p that 

hd t u + \{-hd x + d x V) ■ (hd x + d x V)u = 0, (2.15) 

and here we recognize the Witten operator W = (—hd x + d x V) ■ (hd x + d x V) 
of the introduction. The supersymmetric structure of W is given by 

A = ^l n , <p(x)=i/,(x) = V(x). 



The Kramers- Fokker-Planck operator. The stochastic differential equa- 
tion associated to this model is 

dx = ydt < 21Q 
dy = —'yydt — d x V(x)dt + y/^yhdw. 

The parameter 7 > is a friction coefficient, and the particle of position 
x G R n and velocity y e R n is submitted to an external force field derived 
from a potential V, with w being an n-dimensional Wiener process of mean 
and variance 1. The corresponding equation for the particle density (I2.13P 
is then 

hdj - ihdy ■ (hdy + 2y)f + y ■ hd X f - d X V ■ hdy] = 0. (2.17) 

We have P*(l) = 0, and posing <p(x,y) = y 2 /2 + V(x) we also see that 
P(e~ 2ip / h ) = 0. An application of Theorem 11.21 gives that P is semiclassi- 
cally supersymmetric on R 2n . The (non- normalized) Maxwellian is given by 
M(x,y) = e - 2 ^ x ^' h . 

If we write / = Ai 1 ^ 2 u, then (12.171) becomes, 

1 

hd t u + —(—hdy + y) ■ (hd y + y)u + y ■ hd x u — d x V ■ hd y u = (2.18) 

Taking 7 = 2, we have thus arrived at the semiclassical Kramers-Fokker- 
Planck operator K, mentioned in the introduction, 

K = y ■ hd x — d x V(x) ■ hd y — h 2 A y + y 2 — hn. 

The supersymmetric structure of K is given by 

A = \ ( _i n ^ ) and V&V) = y 2 l^+V(x) 1 *l>{x,y) = y 2 /2+V(x). 
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Chains of oscillators. This is a model for a system of particles described 
by their respective position and velocity (xj, jjj) G R 2n corresponding to two 
oscillators. We suppose that for each oscillator j G {1,2}, the particles are 
submitted to an external force derived from an effective potential Wj(xj), 
and that there is a coupling between the two oscillators derived from an 
effective potential 5W(xx,x 2 ). We denote by W the sum 

W(x) = W 1 (x 1 ) + W 2 (x 2 ) + SW(x u x 2 ), 

where x = (x\, x 2 ), and we also write y = (yi, y 2 ). By Zj, j G {1, 2} we shall 
denote the variables describing the state of the particles in each of the heat 
baths, and set z = (z\,z 2 ). We suppose that the particles in each bath are 
submitted to a coupling with the nearest oscillator, a friction force and a 
thermal diffusion at a temperature Tj = ctjh/2, j = 1, 2. We denote by Wj, 
j G {1,2}, two n-dimensional Brownian motions of mean and variance 1, 
and set w = (wi,w 2 ). The stochastic differential equation for this model is 
the following, 

dx\ = y\dt 

dyi = —d xl W(x)dt + (z± — x\)dt 
I dz\ = —^Z\dt + jxidt — yj^aihdwi , , 

j dz 2 = —^Z\dt + ^x 2 dt — ^J^a 2 hdw 2 

dy 2 = —d X2 W(x)dt + (z 2 — x 2 )dt 

dx 2 = y 2 dt. 

The parameter 7 > is the friction coefficient in the baths. The correspond- 
ing semiclassical equation ( 12. 1 3[) for the particle density is then 

x 2 )/a 2 )f (2.20) 

Z) ■ hdy)f = 0. 

We have -P^(l) = 0. In order to exhibit a semiclassical supersymmetric 
structure for P w , we need to find a smooth function (p(x; h) = <fo( x ) + 0{h), 
such that P w (e- 2 ^ h ) = 0, for aU h > small enough. With a general con- 
figuration, corresponding to different temperatures and non-trivial coupling 
5W 7^ 0, the existence of a Maxwellian of the form e~ 2v ^ h , with ip as above 
is not clear, and in fact we shall show in the next section that in some cases, 
there is no such a Maxwellian, thus establishing part ii) of Theorem 11.31 In 
the remainder of this section, we shall be concerned with two cases for the 
model (12 .20 p . where a supersymmetric structure can be found. 



hd t f + P w f := hd t f + ^a x {-hd Zl ){hd Zl + 2( Zl - 
+ l a2 (-hd Z2 )(hd Z2 + 2(z 2 - 
+ (y ■ hdj - (d x W(x) + x - 
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Let us set W (x) = Wi(xi) + W 2 (x 2 ). The smooth function 

M*,V,*) = E ^ (f + + • (2.21) 

will be of central importance in the following discussion. 

Proof of i) Theorem \1.3\ in the case of equal temperatures. Let us discuss 
first the case when a\ = a 2 =: a. In this case it is immediate to check that 
if we define 

2 1 /,,2 ^ _.^2> 



j=l \ / 



then we have Pw{&~ 2ip ^ h ) = 0. An application of Theorem 11.21 then gives 
that the operator Pyy is supersymmetric on R 6n , in the semiclassical sense. 
The associated Maxwellian is then defined up to a constant by 

M(x,y,z) =e~^ /h . 

Let us also give the conjugated version of the time- dependent problem. If we 
write 

f = M 1 *u, 

the equation f)2.20p becomes 

hd t u + -a (-hd Zl + ^(zi - (hd Zl + ^(zi - xi)*j u 

+ -a (-hd Z2 + ^(z 2 - x 2 )*j (hd Z2 + i(z 2 - x 2 )*j u 
+ (y ■ hd x - (d x W(x) + x - z) ■ hd y ) u = 0. 
The conjugated operator, 

2 



P - 7 

Pw ~ 2 



+ y ■ hd x - (d x W(x) +x-z)-hd y (2.23) 

A,* j 



then enjoys a supersymmetric structure of the form Pw = d^'Td^h, with 



V? = + -<W = - + y 2 /2 + (^ - x) 2 /2) , (2.24) 

a a 
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and a non-degenerate (constant) matrix A given by 



a 



In 

= ^ I ~ l n 

o 7 I n 



(2.25) 



Proof of i) Theorem 11.31 m £/ie case w/ien oW = 0. We consider now 
the case when = 0, so that = Wo, while «i and a 2 m &y be unequal. 
In this case it is immediate to check that with </?o defined in ( I2.2ip . we have 
-2ip Q /h^ _ q Following Theorem ll.2[ we conclude therefore that the 



Pw (e~ 

operator Py/ is supersymmetric, in the semiclassical sense. The Maxwellian 
associated to the problem is then defined up to a constant by 



M(x,y,z) = e ~ 2ipo/h . 



If we write as before, 



f = MV\ 
then the equation ( I2.20P becomes 



7 



hd t u + — oli I —hd zi H (z\ — x\ 



a i 



H — a 2 
2 



-/i<9 22 H (z 2 -x 2 , 

a 2 



/icL, 



CCi 



— ^2 - ^2, 

«2 



(2.26) 



+ (y • — (d x W + x — z) ■ hd y ) u — 0. 



We then have a supersymmetric structure for the conjugated operator, oc- 
curring in (12.26 p . with 



v = i, = ^ = Y^— (Wj(x) + y]/2 + ( Zj - x 3 f/2) , 

3=1 ^ 

and a non-degenerate (constant) matrix A given by 



(2.27) 





/ 








Oil 








o \ 















a 2 








1 




— OL\ 

















2 







— a 2 






























7«i 







V 

















7«2 / 
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3 Non-existence of supersymmetric structures 
for chains of oscillators 



The purpose of this section is to discuss the question of non-existence of 
smooth supersymmetric structures for the model of a chain of anharmonic 
oscillators, in the case when the temperatures of the baths are unequal and 
the effective interaction potential SW is non-vanishing. This will establish 
part ii) of Theorem 11.31 

Following (I2.20p . we are interested in the equation, 



(hdt + Pw)f(t,x,y,z) = 0, t>0, (3.i; 



where 



Pw = iy^ a i(~ hd zj) \ hdz J + ~( z i - x i) ) +y-hd x -(d x W(x)+x-z)-hd y . 
j—i \ j / 

(3.2) 

Here the equation P w (l) = is satisfied automatically. In order to establish 
the breakdown of the semiclassical supersymmetry for Pw, we only have to 
show that there does not exist a function cp e C°°(R 6n ;R) with cp(x;h) = 
(f (x) + 0(h) in the C°°-sense, such that P w (e^ /h ) = for all h > small 
enough, for a suitable choice of the effective potentials Wi, W2, and 5W. 

In what follows, when W = Wo, we shall find it convenient to work with the 
conjugated operator 

P = P Wo := e^ h o P Wo o e-^' h , (3.3) 
where we recall that the function tpo is given by 



j= i J 



a 3 \ 2 

Following (I2.26p . we obtain that 

2 

j — ^ \ 3 y \ 3 

+ y ■ hd x - (d x W (x) + x - z) ■ hd y , (3.4) 
with the function e~ v °l h being in the kernel of Pq and its adjoint. 
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The leading part p = po(x, y, z, £, T), C) of the semiclassical symbol of P — 
P (x, y, z, hD x , hD y , hD z ; h) is given by 

Po = | a i + ( x i ~ z if^ + i V-£- i(d x W (x) + x - z) ■ 7). (3.5) 

We may notice here that Rep > 0. Associated to the operator P in (13.4)) 
is the real symbol 

Qo(x, y, z, f , 77, C) = -Po(a?, V, z, if, ir), <), 

with 

P (x, y, z, hD x , hD y , hD z ) = -q (x, y, z, -hd x , -hd y , -hd z ) 
to leading order. We have 

Qo = ~^2 0Lj (^Cj ~ ^Oj - Zj) 2 J + y-£ - (d x W Q (x) +x-z)-rj. (3.6) 

The phase function ip E C°°{Il en ) satisfies the eikonal equation 

q (x, y, z, d x ip , d y ip , d z (p ) = 0, (3.7) 

reflecting the fact that Po(x, y, z, —ihd x , —ihd y , —ihd z ; h) (e~ <p °l h ') = 0. 

Let Xq E R 2n be a non-degenerate critical point of W (x), and let p E T*R 6n 
be the corresponding point in the phase space, given by x = Xq, y = 0, z = Xq, 
£ = V = C = 0- Let furthermore F Po = F Po , p be the fundamental matrix of the 
quadratic approximation of po at the doubly characteristic point p. Then, as 
explained in [10] , [11] , the spectrum of F po avoids the real axis and is of the 
form ±Afc, 1 < k < 6n, ImA^ > 0. 

Remark. The preceding result comes from the fact that near p, the average 
of Repo along the Hamilton flow of Impo has a positive definite quadratic 
part in its Taylor expansion at p, see [10], [13]. Following [13] and [19], this 
can also be interpreted by saying that the so-called singular space S, defined 
by 

00 

S = ( p| Ker[Re F Po (Im F Po Y]) f|T*R 6 ", 

j=0 

is trivial in this case. 
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Let next p G T*R n and F po be as above, and let F qo be the fundamental 
matrix of the quadratic approximation of go & t the point p. Following e.g. 
|10j . we see that the eigenvalues of F qo and i~ 1 F Po are the same and are of 
the form ±4Afc, 1 < k < 6n, where 

Re [j^kj > °> 1 <k<6n. 

As explained in [10], an application of the stable manifold theorem allows us 
to conclude that in a suitable neighborhood of the point (xo,0,xo) = 7r(p), 
the eikonal equation 

q (x, y, z, d x ip, d y ip, d z <p) = 0, (3.8) 

has a unique smooth solution tp(x,y,z) such that (p{ji(p)) = ip'(n(p)) = 
and (p"(7r(p)) is positive definite. Here n((x,y, z,^,rj,Q) = (x,y,z). Indeed, 
the function tp(x, y, z) is obtained as the generating function for the unstable 
manifold through p for the H qo -Row. Applying this discussion to the case 
when xq is a non-degenerate local minimum of Wq(x), we see that necessarily, 
we have the equality <p = ifo in a neighborhood of 7r(p). 

We shall assume throughout this section that Wo(x) = Wi(xi) + ^2(^2), 
where W2 is a positive definite quadratic form on R". As for W\, we assume 
that Wi is a Morse function on R™ such that d^W x = 0(1) for all a G N n 
with |«| > 2, and such that |VWi(xi)| > 1/C for \xx\ > C > 0. Assume 
furthermore that W\ is a double well potential, so that it has precisely three 
critical points: two local minima m± and a saddle point Sq of signature 
(n — 1, 1). The critical points of Wq are then the local minima M± = (m±, 0) 
and the saddle point So = (so,0), of signature (2n — 1,1). Furthermore, 
W (x) — > +00 as \x\ — > 00. From f)2.2ip we see that the restriction of <p to 
the subspace 

L = {(x,y,z)eR 6n ; z = x, y = 0} 

can be identified with Wq and that the Hessian of </?o in the directions or- 
thogonal to L is positive definite. Consequently, <po is also a Morse function 
on R 6n , tending to +00 when (x,y,z) — > 00, with precisely three critical 
points, given by the local minima M± = (M±, 0, M±) and the saddle point 
So = (So, 0, So), of signature (6n — 1,1). 

Let 5W(x) G C°°(R 2n ), and thinking about SW{x) perturbation, let us 
introduce the perturbed effective potential W = Wo + 5W. The correspond- 
ing operator P\y in (13. 2p is then of the form, 

Pw = Pw - d x (5W(x)) ■ hd y , (3.9) 
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and after a conjugation, we obtain the operator 

P w ■= e ^ /h oP w o e~^l h = P - e w/h o d x (6W(x)) ■ hd y o e~^l h 

= Po-d x (5W(x))-(hd y -d y tpo). (3.10) 

Associated to the operator (13.101) . we have the perturbed real Hamiltonian, 

q(x,y,z,^,r],C) = q (x,y,z,t;,r),() -d x (5W(x)) ■ (rj + dytpo). (3.11) 

We are interested in the question whether the perturbed conjugated operator 
Pw still possesses a smooth supersymmetric structure on R 6n , in the semi- 
classical sense. According to Corollary 2.5, a necessary condition for that is 
the existence of a smooth solution (p of the eikonal equation, 

q(x, y, z, d x ip, d y ip, d z ip) = 0. (3. 12) 

The following is the main result of this section. 

Theorem 3.1 With a suitable choice of W\ and W2 as above, there exists 
5W e C°°(R n ) with M + and Sq ^ supp(SW) such that the eikonal equation 
f)3.12p does not have any solution ip G C°°(Q), with <p(M + ) = cp'(M + ) = 0, 
(p"(M + ) > 0, for any open set C R 6n such that ft + d f2. Here f2 + is the 
connected component of the set ip^ *((— 00, tpo(So))), which contains M + . 

Corollary 3.2 The perturbed operator Pw in (13. 9p does not possess any 
smooth supersymmetric structure in the open set Q, in the semiclassical 
sense. 

Remark. The idea behind the result of Theorem 13.11 is that the solvability 
of the eikonal equation f)3.12p does not seem to be an issue near the local 
minimum M + , but if we start by solving the problem near this point and try 
to extend the solution, we may run into some trouble when approaching the 
saddle point So. 

When establishing Theorem 13. 1[ it will be convenient to write ip = ipo+ip, and 
to consider the corresponding eikonal equation for ip. A direct computation 
using ([MID , (]MD , and fl3TTTD shows that ip e C°°(fi) should satisfy the 
following equation, 

2 

^ + \ a i^f ~ d * 5W ■ W = 29 * SW ■ d v^- ( 3 - 13 ) 
3=1 
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Here the real vector field v is given by 

v = j(z - x) ■ d z + y ■ d x - (d x W (x) + x - z) ■ d y , (3.14) 

and we notice that v can be identified with the Hamilton vector field H qo 
along A^, via the projection 7r\A VQ - Here the Lagrangian manifold A Vo C 
T*R 6n is given by (£,r),() = (d x tp ,d y tp ,d z tp Q ). 

We shall prepare for the analysis of the eikonal equation (13. 13)) by studying 
the behavior of the integral curves of u, in particular close to the stationary 
points of v. Here we may also notice that the flow of v is complete. 

Without using the assumption that Wo is a direct sum, we see that the vector 
field v vanishes at a point (x, y, z) precisely when 

y = j z = Xy d x W o (x)=0. (3.15) 

Recall that for simplicity we assumed 7 = 1, and let xq be a critical point of 
Wq. Then at the corresponding point (xo,yo,Zo) = (xq,0,xq), the lineariza- 
tion of v has the block matrix 



N = -W"(x ) - 1 1 . (3.16) 




Let A be an eigenvalue of N and let (x y z) 1 be a, corresponding eigenvector. 
Then we get the system 

— Aa; + y =0, 
(W"(x ) + l)x + \y - z =0, 
-x + (l-X)z =0. 

Noticing that A = 1 is not an eigenvalue, we obtain that 

y = Ax, z — (1 — A) _1 x, 

where x satisfies 

W"(x )x = - A 2 ^ x. 

Thus, each eigenvalue w of W"{xq) gives rise to three eigenvalues A of N, 
given by the equation 

F(X) := ^ _ A 2 = w . (3.17) 
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We have 



F'(X) 



G(X) 



(3.18) 



(1-A) 



2 ' 



where G(X) = 1 - 2A(1 - A) 2 . It is easy to see that G(X) > for A < 1, 
and therefore, F(X) is strictly increasing from — oo to +oo on the interval 
(— oo, 1), with -F(O) = 0. Next, G(X) is strictly decreasing from 1 to — oo 
on the interval (1, oo), and therefore, on this interval, the function F has a 
unique critical point m > 1, with F(m) < 0. (Numerical computations show 
that m = 1,5652.) 

The equation ( 13 . 1 7|) can be written as a polynomial equation A 3 — A 2 + (1 + 
w)X — w = 0, so that the three roots A-,- satisfy 



We get therefore the following information about the solutions of 03.17p . 

• If w < F(m), then the three roots are real, and two are > 1. The 
remaining root is negative. 

• If w > F(m), then we have precisely one real root Ai and it belongs 
to the interval (— oo, 1). The other two roots are of the form a ± ib, 
with a, b G R, and from the relation (13.191) we see that 2a + Xi = 1, 
implying that a > 0. 

The main conclusion concerning the three roots Xj is as follows: 

• If w > then all three roots have positive real parts. 

• If w = then one root vanishes and the other two have positive real 
parts. 

• If w < then one root is negative and the other two have positive real 
parts. 

• For each eigenvalue w, the x-component of the eigenvector (x,y, zf of 
N is the corresponding eigenvector of W"(xq). 

As the next step in the analysis of the //-flow, we shall show that (po is strictly 
increasing along the integral curves of v in the region where Wq{x) ^ 0. When 
doing so, it will be convenient to write 



3 




(3.19) 



v = vi(w!,d Wl ) +V2(w 2 ,d W2 ), Wj = (xj,yj,Zj). 
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Using f )2.2ip and (13.141) . we get, still assuming that 7 = 1, 

2 



V 



2 x 

j= i o 

2 2 

2 (fo) = Yl ~^~.( Zj ~ x ^ u i( z j ~ x j)> ( 3 - 21 ) 

3=1 J 



2 2 



^o) = Y — i^ z i ~ x ^)) 2 + & - • v2 M ~ x i)) > ( 3 - 22 ) 



3=1 ■ 



^Oo) = ( —(Vj{zj - x j))^A z j - x i) + — Oj - x j) u j( z i ~ x i) 

,=1 V 



(3.23) 



+— - x i)) v A z i ~ x i) + — (*j - ^'KOi - &i) J • ( 3 - 24 ) 
Here 

- x ) = j(z -x)-y, 
u 2 (z - x) = (7 2 - l)(z -x)-jy + Wq(x). 
Still working in the region where Wq(x) ^ 0, we see that 

1) v(<po) > with the equality precisely when z — x. 

2) If u((p ) = and y ^ 0, then u 3 {<^ ) > 0. 

3) If u((p ) = and y — 0, then i/ 2 (z - a;) 7^ and z/ 5 ((/? ) > 0. 

This leads to the statement that for every compact set where Wq(x) 7^ 
there exists a constant C > such that, 

K<A)) > ^ or z/% ) > 1 or ^ 5 (y9 ) > ^. (3.25) 

Proposition 3.3 For every compact set K C R^ n yz where Wq(x) 7^ i/iere 
exists a constant C > snc/i £/iat 

V? o exptz/(x) - ^o(x) > t 5 /C, x G K, < i < 1/C. (3.26) 
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Proof. Consider the function f(s) = i>(ip )oexp sv(x) along an integral curve 
exptv(x) of v. If f(0) > Const. > 0, we get f* f(s)ds > t/C 
for < t <C 1 and (13. 26p follows (for this value of x). 

If < /(0) < 1 and /"(0) > Const. > 0, then f(s) is a strictly convex 
and non-negative function so there exists a unique point s close to 0, where 
f(t) attains its minimum and we have f(s) > (s — s ) 2 /C. Integrating this 
inequality, we see that J f(s)ds > t 3 /C and (ET2Rj) follows. 

It remains to treat the case when < /(0) <C 1, f"(0) < e <C 1 and 
/( 4 )(0) > Const. > 0. The function g(s) = |(/(s) + /(-s)) is even and has 
the Taylor expansion 

= /(0) + ^"(0) 5 2 + ^/ (4) (0)* 4 + 0( S 6 ), 

and is therefore a smooth and strictly convex function of s 2 , non-negative 
for s 2 > 0. Denote this function by k(s 2 ) and restrict the attention to the 
interval [0,t 2 ]. If k'(0) > 0, the strict convexity implies that k(r) > t 2 /C on 
[0,t 2 ]. If k'(t 2 ) < we have k(r) > (r — t 2 ) 2 /C. In the remaining case when 
fc'(0) < and k'{t 2 ) > there exists r G [0, t 2 ] such that fc(r) > (r - r ) 2 /C 
on [0, t 2 } and this is finally the conclusion in all three cases. 

Thus, g(s) > ±{s 2 - t ) 2 for -t < s < t and by an easy calculation, we 
obtain that 

j f{s)ds = j g{s)ds > i J (s 2 - r ) 2 ds >L, < t < 1. 

In other words, (<^o ° exptz/ — ipo o exp (— ti/))(rr) has a lower bound as in 
(j3.26|) . and we are allowed to vary the point x in a small neighborhood, so 
we get the same conclusion for (ip o exp2t^ — ipo)(x) and after replacing 2t 
by t, for (<p o exptv — <p )(x). □ 

From the statement of Theorem 13. 1[ let us recall the set Q + C R 6n , defined 
as the connected component of the set 

^((-oQjVoCSo))), 
which contains M + . It follows from the property 1) prior to (13.251) that 

exp(»(fi + ) C n + , t<0, 

and Proposition ^. 31 gives that exp (tv)(K) converges to when t — > — oo, 

for every K <e 
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The stable manifold theorem tells us that in a suitable neighborhood of So 
there is a unique curve T (manifold of dimension 1), which is stable under 
the y-flow in the sense that if x G T then exp (tu)(x) converges to So, expo- 
nentially fast, when t — > +00. The set T + := T fl Q + is also invariant under 
the forward z/-flow and is the image of a (connected) curve. For wo G r + , 
let us put j(t) = exp (tv)(wo), t G R. Then 7 is a smooth curve in fl + and 
we have, 

'm+. t -> 

(3.27) 

t -> +00. 

The trajectory 7(i) will play a crucial role in the proof of Theorem 13. 11 

Let us now resume the analysis of the eikonal equation (13.131) . The pertur- 
bation 5W(x) will be chosen so that M + supp (SW), and since we are 
interested in smooth solutions cp of (13.121) . for which ip(M+) = <p'(M + ) = 0, 
<p"(M + ) > 0, as we saw above, we have ip = ipo in a neighborhood of 
M + = (M_|_,0,M + ). We shall therefore study the solvability of the prob- 
lem (I3.13p . assuming that ip = in a neighborhood of (M + , 0, M + ). 

Proposition 3.4 Let SW G C°°(R 2n ) be such that M + £ supp (SW) and 
assume that 5W(xi, X2) is a homogeneous polynomial of degree m > 3 in x 2 - 
Ifipe C°°(fi), for £1+ (I n, satisfies fl3~T31 with ip = near (M+,0, M+), 
then we have in Q +! 

ip(x,V,z) = 0((x 2 ,y2,z 2 ) 2 ). 
Proof. We shall view (13.131) as a Hamilton- Jacobi equation of the form, 

where 

P(x, y, z, £, rj, C) = y ■ £ + 7(z - x) ■ ( - (d x W (x) + x - z) ■ rj 



+ \H a ti - dJW ' ^ - 29 JW ' ^0. (3.28) 



In what follows, it will be convenient to write w = (wi, w 2 ), Wj = (xj, yj, Zj), 
and u = (cui, cu 2 ), ujj = (£j, r]j, Q), for j = 1,2. We know that the Lagrangian 
manifold = {(w , ip' (w))} is the if p -flowout of the set 



A^, n neigh((M+, 0, M+; 0, 0, 0),T*R' 
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To be precise, let p(0) = 0(0); w(0)) G neigh((M+, 0, M+; 0, 0, 0), T*R ) be 
such that cj(0) = 0, so that p(0) G C p^ 1 (0), and consider the corre- 
sponding Hamiltonian trajectory 

(w(t); u(t)) = p(t) = exp (tHp)(p ) € A^. 

We shall be interested in the trajectories p(t) for which w 2 (0) = 0. It follows 
from flZZH) and fl3~28j) that 

d W2 p(w, u) = 0{{w 2 , w 2 )), d^piw, u) = 0((w 2 , u 2 )), 

and hence the Hamilton equations 

«>2(*) = ^2^(^(0)^(0)) w 2 (t) = -9 W2 p(w(t),a;(t)), 

imply that 

(«) 2 (t) 7 cJ 2 (t)) = 0(M*W*))). 

Along a trajectory, for which 102(0) = u 2 (0) = 0, we have therefore w 2 (t) = 0, 
cu 2 (t) = 0. A straightforward computation shows next that along the set 
where w 2 = 0, u 2 = 0, we have 

d m p = 0(ui), 

while 

d^p ■ d Wl = viiwi, d m ) + 0(ui) ■ d Wl , 

where u^Wx, d Wl ) = 7O1 - x x ) ■ d Zl + y 1 ■ d xi - [d Xl W x (xi) + %\ - z x ) ■ d yi . 
It follows that an i? p -trajectory pit) = (wi(t),w 2 (t);uJi(t),u 2 (t)) for which 
w 2 (0) = 0, u(0) = 0, satisfies 

pit) = (exp (^(^(O^O^O). 

From the definition of we know that 

dip = udw, (3.29) 

so that ip\w 1 (t),w 2 (t)) = (wi(t),w 2 (*)) = (0,0). Thus, ip'iw^O) = for all 
W\ such that (wi, 0) G Now from (13.291) . by a classical formula given, for 
instance, in Chapter 1 of [5], we also have 

i/>(w(t)) = i/>{w(0)) + [ w(s) • d u p(w(s),u)(s)) ds, (3.30) 

Jo 

where we know that ip(w(0)) = 0, so that tp(wi, 0) = again for all w\ such 
that (u>i,0) G Vt + . Using this and ip'(wi,0) = 0, we get that ip(w) = 0{w 2 ) 
and the proof is complete. □ 
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In what follows, we shall assume that SW G C°°(R n ) has the properties 
described in Proposition 13.41 Using also that 

d yj <p = —yj, i = 1,2, 

we see that the right hand side of ( 13 . 13[) is homogeneous of degree m in 
W2 = {x-iiViiZ-}). If the problem (13. 13[) has a smooth solution ip G C°°(Q), 
then we have a Taylor expansion at W2 = 0, writing also w\ = (xi, 2/1,21), 

CO 

ij){w u w 2 ) ^ ^{wi, w 2 ). (3.31) 

k=0 

Here ^(wi, u>2) is homogeneous of degree k in w 2 , with C°°-coemcients in 
w\, and by Proposition I3.4[ we know that ipo and ipi vanish in Q + . It follows 
from f)3.14p that uipk is also homogeneous of degree k in W2- We see then 
that the term, homogeneous of degree fi > 0, in the left hand side of (I3.13p . 
is given by 

fc=0 k=0 

- d xl 5W ■ d yi tp d X2 SW ■ (3.32) 
Here it is understood that ipj = 0, for j < 0. 

It is now easy to conclude that ip^ all vanish, for /i < m, on the open set 
f2 + C R 6n . Indeed, taking [x = 2 in (I3.32p . we see that the sum 

2 

k=0 

vanishes in Q + , while the only non- vanishing term in the sum 

2 

^2(d Z2 ifj k+1 )(d Z2 ij2-k+i) 

k=0 

is given by (d Z2 ip2) 2 - We get the equation 

Hi + \u2(d Z2 i)2f = 0, (3.33) 

where we also know that ip 2 vanishes near (M+, 0, M + ). The preceding equa- 
tion can be viewed as a first order ordinary differential equation for ip 2 along 



28 



the integral curves of v in Q + n {(u>i,u> 2 ) ; w 2 = 0} and we can conclude 
that ^ 2 = in 

We shall now see, arguing inductively, that ijj^ — in for /i < m. Indeed, 
assume that ip = ipi = ip2 = • • • = VV-i = 0' f° r some 2 < /i < m. Then we 
have 

k=0 

As for the sum 

fc=0 

we see that the only term here which is not clearly vanishing corresponds to 
the case when k + 1 = fi. In this case, the corresponding term is equal to 
(d Z2 i/j^)d Z2 ip2, which vanishes after all. The sum above consequently vanishes, 
and from (I3.32p we get the equation 

vipfj, = 0, fi < m. 

It follows that ipn = in Q + , for fi < m. We conclude that a smooth solution 
ip of (I3.13p . such that ip = near (M + , 0, M + ), has the following form in Q + , 

^(w u w 2 )=^ m (w 1 ,w 2 ) + 0(w 2 n+l ). 

Combining (13.131) and (13. 32p . we see that ip m should satisfy the following 
non- homogeneous transport equation in Q + , 

vi) m = 2dJW ■ d y ipo. (3.34) 



The proof of Theorem 13.11 will therefore be concluded, once we establish the 
following result. 

Proposition 3.5 There exist a positive definite quadratic form W^O*^)? a 
Morse function Wi(xi) with two local minima and a saddle point, and a 
perturbation 5W e C°°(K 2n ) with M + supp (SW), S £ supp((W), and 
such that 5W(xi,X2) is a homogeneous polynomial of degree m > 3 in x 2 , 
for which the transport equation (I3.34p does not have a smooth solution in 

n + u{s }. 

Proof. With the notation Wj = (xj,yj,Zj), j = 1,2, let us write v = 
Y^j=i u j( w jt9 Wj ). The preceding discussion shows that there exists an inte- 
gral curve 71 of V\ such that 

j(m+,0,m + ) t -00, 
l(s ,0, s ) t y +00. 
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Let N 2 be the coefficient matrix of the linear vector field v 2 , which we shall 
view as the linearization at w 2 = 0, and let Ai, A3„ be the corresponding 
eigenvalues, so that Re Xj > 0, 1 < j < 3n. Let us assume, as we may, 
that N 2 has no Jordan blocks, and after a linear change of variables, we may 
therefore assume that w 2 = (cui, ...,cjs n ), and 

3n 

Then, writing A = (Ai, . . . A 3n ) G C 3n , we have 



Let us consider the equation (I3.34p . 

(yx + Hj)(Vm) = —Vx ■ d Xl SW + — y 2 ■ d X2 SW, (3.35) 

and put 

2 

ip m = — SW(x) + u. 
a.\ 

Then we get 

{ Vl + u 2 )(u) =(^-l-)y 2 . d X2 5W, (3.36) 
\a 2 ax) 

where we recall that ot 2 7^ a\. Using that 5W(xi,x 2 ) is a homogeneous 
polynomial of degree m > 3 in x 2 for every x\, we may write 

( — - — ) V2 ■ d X2 5W = V g a (x!)uj a , 



\a 2 a\ ) 



\a\=m 



and if (I3.36P has a smooth solution u, we can assume without loss of generality 
that 

u= } j u a (wi)u a . 

\a\=m 

The equation H3.36|) reduces then to the following decoupled system of equa- 
tions, 

{yi{w u d Wx ) + A • a)u a (wi) = g a {xx), \a\ = m. (3.37) 

We shall choose 5W(x 1 ,x 2 ) = ±ip(xi)v(x 2 ), where < ip G C^R"), 
m + ,SQ supp^, while ip > somewhere on the image of 71, and where 
v is a homogeneous polynomial of degree m. Then for some a = ao of length 
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m, we have < g a G C^°(R rt ), m + ,s $ supp^, and g a is > somewhere 
on the image of 71 . 

The equation (I3.37P along jx with a = a has a compactly supported right 
hand side with constant sign, not identically equal to 0, so the solution 
u — u ao is non-zero either on 71 D neigh (m + , , m + ) or on 71 D neigh (s , , Sq) . 

In the first case we use the natural parametrization 71 (t) = expt^i(wo), 
—00 < t < +00, where wq is some fixed point on 71. Then there is a 
constant C > such that dist (71 (t), (m + ,0,m + )) < Ce~^ /C for t < 0. For 
t large and negative, we have 



Thus u a is unbounded near (m + ,0,m + ), and hence cannot be smooth near 
that point. 

In the second case, we recall that 71 is a part of the one-dimensional stable 
manifold through (so, 0, sq) for the 1/1-flow. We can find new smooth local 
coordinates x = (xi,x") centered at (sq, 0, sq), such that this stable manifold 
is given by x" = 0, and hence 



where aj(xi,0) = when j > 2. Furthermore, ai(xi,0) = — jii(xi + f(xi)) 
where fix > 0, f(xx) = 0(xf) and we may assume that 71 coincides with the 
positive Xi-axis near (so,0, sq). Along 71 and near (so,0, s ) we know that 
u a is a non- vanishing solution of the following equation, 




and hence, 



3n 




(-//i^a + f(xx))d xi + A • = 0, < xx < 1, 



so that 




) 



where C ^ and > is small and fixed. Here 



1 



1 1 



1 m , m 2 



s + f(s) 




s s 2 s 3 
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so 

(•XI 



f X1 1 

/ — , , ds = In xi + g(x 1 ) 



z° S + f{s) 

where g is smooth near x\ = 0. Thus, 

u a (xi,0) = Cx{?e 9ixi \ C ^ 0, 

and if ^ ^ N (which can be arranged by choosing the parameters suitably, 
c.f. ( I3.17P ). we conclude that u a cannot be smooth near (s , 0, s ). The proof 
of Proposition E3] is complete. □ 
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